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' Abstract. The presence of electron-electron interactions in one dimension profoundly changes 

the properties of a system. The separation of charge and spin degrees of freedom is just one 
example. We consider what happens when a system consisting of a ferromagnetic region of 
1) ' non-collinearity, i.e. a domain wall, is coupled to interacting electrons in one-dimension (more 

»' | , specifically a Luttinger liquid). The ferromagnetism breaks spin charge separation and the 

presence of the domain wall introduces a spin dependent scatterer into the problem. The 
absence of spin charge separation and the effects of the electron correlations results in very 
. different behaviour for the excitations in the system and for spin-transfer-torque effects in this 

model. 
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i ^i ' 1. Introduction 

The behaviour of the electronic and magnetic degrees of freedom in quasi one-dimensional wires 
and films has already received considerable interest and study, see for example the review of 
■ Marrows[T]. So far most works focus on how the transport properties of free electrons behave in 

a ferromagnetic wire with a domain wall, and how these spin polarized currents set the domain 
wall itself into motion. One obvious question to ask is whether interactions are important in 
such cases. A straightforward yes was answered by Dugaev et al.[2] who explicitly demonstrated 
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build up around the domain wall. This of course should not be surprising but it is a timely 
reminder that correlation effects are often important even for nominally "free" electron systems. 
In this case failure to include these correlation effects results in an un-physical charge excess. 
In addition to the mean field interaction work of Dugaev et al. there also exists in the literature 
some consideration of Hartree-Fock corrections [SJ 3]. Here we will consider a more strongly 
correlated system: the Tomonaga-Luttinger liquid. 

Experimentally [S] the construction of chains of single magnetic atoms is already possible and 
one dimensional magnetic systems are experimentally available. Several obvious questions must 
be answered at this point. Traditionally it is thought that long range ferromagnetic order does 
not exist in one dimension. However, there are several things which mitigate this bald statement 
for these systems. The long range ferromagnetic order is in fact broken by the presence of domain 
walls, precisely the situation we are interested in. We also note that effects from the substrate 
mean that one can not really consider the chain as an isolated spin chain to be solved by a 
Heisenberg or Ising type model. Experimentally ferromagnetically ordered chains of atoms are 
indeed observed [5j [61 [8]. Furthermore they exhibit non-collinear ferromagnetic order[9j. One 
should of course bear in mind that the condition for defining the system as one-dimensional may 
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be different for the magnetic degrees of freedom and for the electronic degrees of freedom. It is 
only the latter that we require to be one-dimensional to be in the Luttinger regime. 

These considerations lead us to our model of a ferromagnetic Tomonaga-Luttinger model[10, 
[TTJ [TJl [13] . Unfortunately the chains of atoms do not easily display Luttinger physics due to the 
interference of the substrate [2], though they remain an intriguing possible future application 
of the model in a more complicated scenario. Nonetheless there is still a class of possible 
candidates left for our model: dilute magnetic semiconductors [J5]. Systems where the magnetic 
and electronic degrees of freedom belong to different layers would also be a possible realization. 

We will consider an "s-d" like model in which the bulk magnetization and the conduction 
electrons are treated as separate (though of course still interacting). As such we introduce 
two timescales into the problem, a fast electronic one and a slow magnetic one. This allows 
us to first answer the question of how the presence of the domain wall affects the Luttinger 
liquid, forgetting for the moment the effect the moving domain wall will have on the conduction 
electrons. 

The question of a non-collinear ferromagnetic Luttinger liquid has not yet been fully 
addressed [13J, we set up the model in section (2) Firstly we would like to analyze the low energy 
effective field theories and their low temperature properties, the results of this are presented 
in section [3j Secondly we wish to understand the behaviour and nature of the appropriate 
excitations of the model. Finally we want to consider how this affects the magnetization 
dynamics of the domain wall. Due to the radically different nature of the excitations of a 
Luttinger, as opposed to Fermi, liquid we expect these dynamics to be very different. We offer 
some thoughts on this based on our preliminary results in section H] later in this article. 



2. The Model 

The magnetization of the wire is described by M(x) = M s cos[0(x)]z + M s sin[0(x)]y where M s 
is the saturation magnetization and, using cos[G(x)] = — tanh[x/A], we have a description of 
a domain wall of length A situated at x = 0. This magnetization is coupled to the conduction 
electrons with a strength given by the exchange coupling J. Naturally we consider a screened, 
and hence short range, interaction V(x-x'). This allows us to start from the following standard 
"s-d" Hamiltonianll6l: 



ipajx) - J y~] J dx^l(x)ijj a f(x)a aa >.M(x) 

era' 

+ J dxdx'ipl(x)i'l,(x')V(x - x')ip a /(x')ijj a (x) . 
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li is the chemical potential. 

In order to be able to linearize the system our first step must be to remove the spatially 
dependent, and in principle perhaps very large, magnetization. To this end we rotate the 
spin direction to get a collinear ferromagnet via the following gauge transformation [T71 I18| : 

H = U t (x)HU(a;), ip a (x) = Ul rr ,{x)^ al (x), and U(s) = eW "* . The interaction is left 
unaffected, the magnetization is locally rotated to a Zeeman term, and the kinetic energy 
operator introduces a new potential when it acts on the local rotation. Thus we have 
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d 2 x - JM s a z aa - ll 



-0o-(x) + Hj + H w and 



d^t(x)[ie'(x)a^d x + l ^o-l a , - \[6'(x)]H aa ,]^(x). 
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This new potential can be approximated if we assume that the Fermi wavelength Af < A. In 
this case only the first term of equation ([3]) is relevant, the next two terms are of order (Xp/X) 2 . 

The next step is linearization via the ansatz if) a (x) = e lkFcrXr ipR a {x) + e~ lkpiyX ^L a (x), where 
kpa = y / 2m(/i + JM s o z aa ). The R and L indices denote right and left moving electrons 
respectively. Note that if the Zeeman term is large we must linearize around the spin split 
Fermi points. This leads to the physically relevant phenomena of the breakdown of spin-charge 
separation [TS]. Umklapp processes scattering two left movers into right movers and vice versa 
are always neglected here due to the non-commensurate nature of the Fermi wavevectors. 

Finally we have our model to be bosonizedQ-OJ. Naturally the bosonic fields satisfy 
[4> (T (x),Il a (x')] = iS(x — x') where H a (x) = d x 9 a (x). We first rotate to the spin and charge 
degrees of freedom. We can then diagonalize the quadratic part of the Hamiltonian with 
(j) c i s (x) = [0i (x) ± <f)2{x)]/y/2 (and similar for the 6{x) fields) giving us for the quadratic part 
of the bosonic Hamiltonian 



H n = / dx 



2V C K C , . 2 v sK.? 



[d x e c {x)Y^ + [d x e s {x)Y^ + [d x e c (x)d x o,(x)]v a 

Md x Mx)] 2 ^ + [d x Mx)] 2 ^ + [d x Mx)d x Mx)]v b • (4) 

K s and K c are the spin and charge Luttinger parameters. v a and v b describe the coupling between 
the spin and charge sectors. These parameters are functions of the interaction strengths and 
Fermi velocities. 

The diagonalization of the quadratic terms introduces new velocities, Ui, for Hq and a set 
of {Tf'^} parameters. These are all known in terms of the previously mentioned Luttinger 
parameters. The final Hamiltonian is H = Hq + Hi + Hw + H^ where 

H = f dx Y, y i(9Ji(x)) 2 + (H 4 (x)) 2 ] (5) 

J i=l,2 

H l = ^J Mvpt + v n\ sed f /A] M^Pf h{x) + Tium 

sin[a;(fc n - k n ) + v^T^x) + T*fa{x)] (6) 

< = w~ I Mm - v Fi ] secI f /A] MV^iT^ix) + Tfe 2 {x))} 
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sm[x(^ t + k n ) - V2^(f?Mx) + TtMx))] (7) 
Hi = j^-^ J dxcos[2v / 2^(T 1 ^ 1 (x)+T 2 ^ 2 (x))] 

[[ffiJ-t + 9i±i\ cos[2x(k Ft - k Fi )} - % 1±t - g 1±i ] sin[2x(A; Ft - k Fi )}] . (8) 

We have both forward and backward scattering terms, Hi, and H^ respectively. We also find 
an oscillating sine-Gordon interaction term Hi, describing "slow" Umklapp processes which are 
not averaged out, with some strength gi± a . These all couple our otherwise diagonal bosonic 
degrees of freedom. The appropriate excitations of such an SU(2) asymmetric model have no 
obvious physical interpretation. This summarizes the effective bosonic field theory which is the 
foundation of our mathematical analysis of the problem. 

3. Low Energy Physics 

We start by writing a functional integral partition function [20J: 



with periodic boundary conditions in imaginary time r. Following the standard procedure we 
split the fields into fast, and slow, ^ < , fields. Our fast fields are defined for A' < |A;|, \u\/ui < 
A, and the slow for |fc|,|u;|/uj < A'. Expanding the exponent in terms of H\ and H w and 
performing the averaging over the fast modes we then re-exponentiate the expression to find the 
appropriate scaling equations. Parameterizing A^'^ = Aoe~ l (~ dl ^ we find for the first order Hi 
term 

~1T = 91± ^ ■ (10) 

7 1 , and also the yf ,b used below, are known functions of the rotation terms {T- '^}. This term is 
an irrelevant perturbation for any situation we wish to look at. In the limit of weak magnetization 
we can simplify the expressions to find 7 1 — > 2 — 2K S . In this limit K s > 1 and it becomes clear 
that the term is irrelevant. 

The second order equation for our model is more difficult than for the sine-Gordon model. 
A diagonal equation in the (pi's is not recovered and to perform any further analysis we would 
have to rediagonalize the problem and then renormalize the model once again. This is perhaps 
not totally unexpected as the scattering term we are dealing with explicitly couples these terms. 
We leave the more involved second order renormalization group analysis to a future work and 
stick here to the first order equations. 

The same analysis is performed on the scattering terms, we define the scattering coupling 
constants g^ b (0) = vf^ ± vp±. For back scattering we find: 

This term is a relevant perturbation. In the limit of weak magnetization this simplifies to 
7 6 — > 2 — [(ET S ) _1 + K c ]/2. Similarly the forward scattering equation stands as 

< 12 > 

This is also a relevant term. In the limit of weak magnetization this becomes 7^ — > 2 — [(ET S ) _1 + 
K s }/2. 

In contrast to the scaling equations of Pereira and Miranda|13t |2~T] we find that we can not 
neglect forward scattering. They were interested in the limit of very sharp walls and the small 
magnetization limit. The model they solve leads to the consideration of a magnetic impurity 
in an otherwise collinear ferromagnetic system. In such a case it is perhaps not surprising 
that this model leads to the magnetic analog of the Kane-Fisher problem |22[ 123] . Their final 
bosonic Hamiltonian is appropriate only for very sharp walls where the domain wall profile will 
be different to the longer walls which we wish to consider, which leads to the different scaling 
equations and hence a different effective low temperature model. It should also be pointed out 
that depending on the length of the domain wall forward scattering can become the dominant 
form of scattering. Of course in the physical model "forward scattering" refers to an electron 
which passes through the wall without changing its spin. 

We have three natural length scales present in the problem: A+ = 2-K{kp\ + kpi)" 1 , 
A_ = 2ir(kFX ~ and A. In the limit of weak magnetization A_ — > 00 and 2A+ — > A_p and 

the limit of large magnetization, when one spin channel becomes frozen out gives A+ ~ A_ ~ Xp. 
We are generally interested in long domain walls when A > Xf and we are in the physical regime 
where both forward and backward scattering are important. As one might expect in the large 
magnetization limit the scattering terms can be neglected as only one spin channel is available. 
For very sharp delta function like walls one finds only back scattering relevant [13] , a case we do 



not consider further. For long domain walls we are back in the adiabatic limit where scattering 
can be completely neglected. Interestingly there is a regime in between, A + < A < A_, where 
forward scattering is the dominant mode of scattering. This regime can also be reached by 
tuning the magnetization to be suitably weak. 

4. Outlook 

We are principally interested in the spin density and the spin and electronic currents. The 
spin density and current will also necessary for understanding the magnetization dynamics of 
the domain wall. The spin density [241 [25] is defined as S(x) = (x)U^ aUij>{x) . It can 
be calculated perturbatively in the domain wall profile using the bosonic model derived above. 
Preliminary results show an important development. The corrections to the zeroth order density 
caused by scattering are concentrated around the edges of the domain wall. In a normal quasi 
one-dimensional wire these corrections are found around the centre of the domain wall [2] and 
contribute to a spin torque which causes precessional motion of the domain wall, in analogy with 
a moment in a magnetic field. In contrast we expect the domain wall to be distorted rather than 
set into precessional motion by these terms in the Luttinger liquid. The temperature dependence 
of the spin density corrections can be extracted from these corrections for the small temperature 
limit. One finds power law contributions in the forward and backward scattering exponents: 
T~"i b . These become sharper for larger values of JM. 

The magnetization dynamics are described by the Landau-Lifschitz-Gilbert equation, or some 
suitable generalization thereof |26j [271 128] . There are two different aspects to this. One is the 
straightforward point that the dynamics, over the length of the wall, will be affected by the 
different spin density of the Luttinger liquid compared to the Fermi liquid or non-interacting 
case. The second, more interesting point, is whether the derivation of the non-adiabatic terms 
in the LLG equation are valid for a Luttinger liquid. 

Following Zhang and Li [28] one can derive contributions to the magnetization dynamics which 
allow for the fact that the electrons do not instantaneously follow the magnetization profile. One 
first writes a continuity equation for the spins, assuming part of it to be always parallel to the 
bulk magnetization and allowing a small deviation from this. In order to derive the current 
dependent (so called (3-) terms, those which drive the domain wall along the wire, one assumes 
that js' y ' z (x,t) = —fj,BPje' y ' z (x,t)M x ' y,z (x,t)/eM s . j e (x,t) is the charge current and P is the 
magnitude of the polarization, whilst j s (x, t) is the spin current. A quite reasonable assumption 
in a Fermi liquid, this of course starts to look more dubious in the case of a Luttinger liquid. 
When SU(2) symmetry is present, not our case one is reminded, spin and charge are of course 
uncorrelated and possess different velocities. Thus this assumption would completely fail. For 
us the situation is not so simple, nonetheless what is obvious from the model is that spin and 
charge are not fully correlated. One is forced to work with the spin current and not the electric 
current and, as we have already seen, the spin degrees of freedom can behave rather differently 
for this model. 

Other than the magnetic behaviour we also wish to analyze the transport through the wire. 

5. Conclusion 

We have considered a non-collinear ferromagnetic Luttinger liquid. We assessed its low 
energy properties and the relevant lengthscales, presenting the results of renormalization group 
calculations. It is possible to tune the relevant types of scattering operators by considering 
different domain wall lengths and magnetization magnitudes. Here 'magnetization' refers to 
both its magnitude and the magnitude of the exchange coupling. The regime of most interest to 
us which we go on to assess is the intermediate regime where we have both forward and backward 
scattering terms and an appreciable magnetization. We discuss some preliminary results for the 
spin density correlation function and the effect that this will have on the subsequent domain 



wall dynamics. We propose that in our Luttinger liquid system the magnetization dynamics will 
be altered by the very different nature of the Luttinger liquid excitations. 
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